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Detection of superhumps in XTE J1118+480 approaching 
quiescence 
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ABSTRACT 

We present the results of our monitoring of the halo black-hole soft X-ray transient 
(SXT) XTE J1118+480 during its decline to quiescence. The system has decayed 0.5 
mags from December 2000 to its present near quiescent level at R~18.65 (June 2001). 
The ellipsoidal lightcurve is distorted by an additional modulation that we interpret 
as a superhump of P s h=0. 17049(1) d i.e. 0.3% longer than the orbital period. This 
implies a disc precession period P pre c ~ 52 d. After correcting the average phase-folded 
light curve for veiling, the amplitude difference between the minima suggests that the 
binary inclination angle lies in the range i = 71 — 82°. However, we urge caution 
in the interpretation of these values because of residual systematic contamination of 
the ellipsoidal lightcurve by the complex form of the superhump modulation. The 
orbital-mean Ha profiles exhibit clear velocity variations with ^500km/s amplitude. 
We interpret this as the first spectroscopic evidence of an eccentric precessing disc. 
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1 INTRODUCTION 

XTE J1118+480 is an important SXT for a number of rea- 
sons. It is at high galactic latitude (6 = 62.3°) which makes 



it the first halo black-hole bi nary, since it mnsfr be~at 



with reprocessing in the accretion disc ( Haswell et al. 2000 ) . 
The optical rise to outburst preceded the X-ray rise by ~ 
10 d which allowed the size of the advective c orona to be 
constrained to ~ 1.2 x 10 4 Schwarzschild radii (Wren et al 
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et al. 2001). The low interstellar absorption allows for de- 



2001). The radial velocity curve of the companion star and 
the large mass function (6.1 ±0.3 Mq) demonstrates that 
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XTE J1118+480 contain s a black hole (McClintock et al 



of an optically thick disc plus synchrotron emission, whereas 
the EUV to X- ray spectrum is rem iniscent of the low (hard) 
state of SXTs (Hynes et al. 2000). Since L x /L opt — 5 com- 
pared to typically 500 for galactic low mass X-ray binaries 
(LMXBs), it has been suggested that it may be an Ac cretion 
Disc Corona (ADC) source seen at high inclination ( |G!arcia 
et al. 200C ). However, the detection of X-ray, UV and optical 
QPOs at 0.08 Hz suggest that the inner disc is directly vis- 
ible. The UV variability lags the X-rays by 1-2 s, consistent 
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2001; Wagner at al. 2001). Finally, the rotational velocity of 
t he compani on star implies a mass ratio of q=0.037±0.007 



( |Orosz 200i| ) 



Optical superhumps were observed during the o utburst, 
changing shape and period ( Ucmura et al. 2000a| ). Super- 
humps are optical modulations first discovered in Dwarf 
Novae superoutbursts with a period a few percent longer 
than the orbital period and a non-sinusoidal shape. Since 
then, this phenomenon has been seen in decaying novae, 
AM CVn systems and novalike variables. Both the ampli- 
tude and shape of the superhump modulation change while 
the outburst declines. The most promising models in ex- 
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plaining the superhump behaviour assume that the accre- 
tion disc expands, due to the action of viscosity, to the 3:1 
resonance radius and the eccentric disc is then forced to pre- 
cess by perturbations from the secondary (e.g. Whitehurst 
& King, 1991). Superhump lightcurves can be explained by 
changes in the disc luminosity associated with the periodic 
deformation of the disc shape (e.g. Simpson & Wood, 1998). 
Since the resonance is only possible for small mass ratios, 
superhumps should also appear in SXTs, as has been con- 
firmed by O'Donoghue & Charles (1996). Furthermore the 
strongest resonance (2:1) requires very extreme mass ratios 
(q<0.025) which might be reached in some SXTs such as 
Jl 118+480 as we will discuss here. 

In this paper we present the detection of superhumps 
in XTE J1118+480 when the system was approaching the 
quiescent state. This discovery provides the first solid evi- 
dence of a precessing disc in a SXT near quiescence. Super- 
humps in dwarf novae were traditionally though t never to 
be present in quiescence but recent observations (Patterson 
et al. 1995) have shown that superhumps can indeed persist 
into quiescence after a superoutburst has ended. The case of 
XTE J1118+480 may be analogous behaviour in the SXTs. 
A preliminary analysis of these data was reported in Casares 
et al. (2001). 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Photometry 

Jl 118+480 was observed for a total of 53 nights in the period 
December 2000- June 2001, during the time when the system 
was approaching quiescence, mainly in the 7?-band with six 
different telescopes: 

• the 0.82 m (IAC80) and the 1 m Optical Ground Sta- 
tion (OGS) telescopes, equipped with identical Thomson 
CCD cameras and using exposure times ranging from 300 
to 1200 s , at Observatorio del Teide (Tenerife). The IAC80 
and OGS were operated simultaneously for two whole nights 
in order to obtain 71-7 colour information. 

• the 2.5 m Nordic Optical Telescope (NOT) and the 
1 m Jacobus Kapteyn Telescope (JKT) at Observatorio del 
Roque de los Muchachos (La Palma) . The data at the NOT, 
which have the best time resolution, were taken with AL- 
FOSC with exposure times of 30 and 60 s. At the JKT we 
used exposure times ranging from 90 to 600 s with the SITe2 
detector. 

• the 1.3 m McGraw-Hill Telescope (MGHT) at the 
MDmQ Observatory (Arizona) and the 1.55 m Kuiper Tele- 
scope of the University of Arizona on Mount Lemmon. Pho- 
tometry in the 7 band was obtained with the Kuiper using 
the 2kx2k CCD camera (2kBigCCD) and in R band with 
the McGraw-Hill Telescope and 2kx2k Echelle CCD cam- 
era. The exposure time for both sets was 300 s and the seeing 
was typically 0.8-1.5 arcsec. 

In all cases we used 2x2 pixel binning to improve the 



t The MDM observatory is a joint facility of Dartmouth College, 
The University of Michigan, Columbia University and the Ohio 
State University. 



readout time. In Table |l] we present a log of the observa- 
tions. 

In total, the target was observed for over 70 orbital cy- 
cles. The individual images were de-biased and flat-fielded in 
the standard way, with the data reduction being performed 
within IRAF^. The instrumental magnitudes were obtained 
using PSF ph otometry with the iraf routine daop hot 
( [gtetson 1987 ) or optimal photometry (Naylor 1998) us- 
ing the Starlink photom package. Differential light curves 
were constructed relative to a nearby comparison star. Light 
curves of this comparison star with respect to others were 
also obtained in order to check for variability from which we 
estimate that our differential photometry is accurate to ~1 
per cent in all observations. 



2.2 Spectroscopy 

We observed J1118+480 using the blue channel CCD spec- 
trograph on the 6.5 m Multiple Mirror Telescope^] (MMT) 
at Arizona, on the nights of 2000 Nov 20, 30 and 2001 
Jan 4. The seeing was typically 1 arcsec and a 1 arcsec wide 
slit was employed. Each 1440 s exposure was bracketed by 
a HeNeAr lamp spectrum which led to a wavelength cali- 
bration accurate to 5-7kms~ rms. Spectra were also ob- 
tained on 2001 Jan 12, with the ISIS red channel of the 
4.2 m William Herschel Telescope (WHT) at Observatorio 
del Roque de los Muchachos (La Palma). The slit width 
was 1.5 arcsec and wavelength calibration was checked with 
respect to night-sky emission lines to be within lOkrns -1 . 
Further details of these observations can be found in Wag- 
ner et al. 2001; A further 29 spectra were taken on 15 and 
16 Apr 2001 with the Steward Observatory 2.3 m Bok Tele- 
scope on Kitt Peak equipped with the B&C Spectrograph. 
A 400 line mm -1 grating and 1.5 arcsec wide entrance slit 
were employed which yielded spectra covering the spectral 
region 4100-7400 A at a spectral resolution of 5.5 A. The 
exposure time was 1400 s in typically 1.4 arcsec seeing and 
photometric conditions. The spectroscopic dataset was com- 
pleted with 49 spectra obtained on 27 and 28 Apr 2001 
on the 4.2 m WHT using the ISIS red arm. A narrow slit 
(0.8 arcsec) aligned with the parallactic angle and with no 
comparison star in the slit, was used in combination with a 
1200 line mm -1 grating. The wavelength calibration derived 
from arc lamps was checked against night sky lines and was 
accurate to lkms -1 (see table ^). 

Standard iraf procedures were used to de-bias the im- 
ages and to remove the small scale CCD sensitivity varia- 
tions. One dimensional sp ectra were ex tracted using the op- 
timal extraction method (Home 1986). Where possible we 
performed flux calibration relative to other stars on the slit. 



• IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under contract with the National 
Science Foundation. 

§ MMT is operated as a joint facility of the Smithsonian Institu- 
tion and the University of Arizona by the Multiple Mirror Tele- 
scope Observatory. 
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Table 1. Log of photometric observations 



Night 


Telescope 


Filter 


Exposure time (s) 


Coverage (cycli 


2000 Dec 13 


IAC80 


R 


600 


0.97 


2000 Dec 14 


IAC80 


R 


600 


1.01 


2000 Dec 28 


IAC80 


R 


600 


1.22 


2001 Jan 09 


IAC80 


R 


600 


0.50 


2001 Jan 11 


IAC80 


R 


600 


0.97 


2001 Jan 19 


IAC80 


R 


600 


0.80 


2001 Jan 27 


IAC80 


R 


600 


0.93 


2001 Jan 31 


IAC80 


R 


600 


1.89 


2001 Feb 01 


IAC80 


R 


600 


1.76 


2001 Feb 02 


I AC 80 


R 


1200 


1.58 


2001 Feb 03 


IAC80 


R 


900 


1.76 


2001 Feb 04 


IAC80 


R 


600 


1.35 


2001 Feb 24 


IAC80 


R 


600 


1.01 


2001 Feb 27 


IAC80 


R 


600 


0.34 


2001 Mar 07 


IAC80 


R 


1200 


1.16 


2001 Mar 08 


IAC80 


R 


600 


1.01 


2001 Mar 09 


IAC80 


R 


600 


0.97 


2001 Mar 10 


IAC80 


R 


600 


1.08 


2001 Mar 11 


IAC80 


R 


600 


1.76 


2001 Mar 18 


Kuipcr 


I 


300 


1.50 


2001 Mar 19 


Kuipcr 


I 


300 


1.69 


2001 Mar 20 


Kuiper 


I 


300 


1.37 


2001 Mar 21 


Kuipcr 


I 


300 


1.69 


2001 Mar 25 


IAC80 


R 


600 


1.87 


2001 Mar 26 


IAC80 


R 


600 


1.64 


2001 Mar 28 


I AC 80 


R 


600 


1.91 


2001 Apr 01 


NOT 


R 


30 


2.02 


2001 Apr 02 


NOT 


R 


60 


1.34 


2001 Apr 04 


IAC80 


R 


600 


1.01 


2001 Apr 05 


OGS 


R 


300 


0.83 


2001 Apr 06 


OGS 


R 


300 


1.68 


2001 Apr 07 


OGS 


I 


300 


1.76 


2001 Apr 08 


OGS 


R 


300 


2.23 


2001 Apr 09 


OGS, IAC80 


R, I 


420, 420 


1.85, 1.75 


2001 Apr 10 


OGS, IAC80 


R, I 


420, 600 


2.00, 1.60 


2001 Apr 11 


IAC80 


R 


600 


1.83 


2001 Apr 12 


IAC80 


R 


600 


1.60 


2001 Apr 13 


IAC80, MGHT 


R, R 


600, 300 


1.90, 1.59 


2001 Apr 14 


MGHT 


R 


300 


1.26 


2001 Apr 15 


MGHT 


R 


300 


0.89 


2001 Apr 16 


MGHT 


R 


300 


1.72 


2001 Apr 17 


MGHT 


R 


300 


1.69 


2001 Apr 20 


JKT 


R 


90-120 


0.50 


2001 Apr 22 


JKT 


R 


90 


1.86 


2001 Apr 23 


JKT 


R 


600 


0.27 


2001 Apr 24 


JKT 


R 


90-600 


1.77 


2001 Apr 25 


JKT 


R 


90-180 


1.21 


2001 Apr 27 


JKT 


R 


180 


1.12 


2001 Apr 28 


JKT 


R 


300 


1.70 


2001 May 29 


IAC80 


R 


600 


0.88 


2001 May 30 


IAC80 


R 


1200 


0.91 


2001 Jun 12 


OGS 


R 


600 


0.64 


2001 Jun 26 


IAC80 


R 


600 


0.25 



\es) 



Absolute photometric calibration was not attempted due to 
slit losses. 

3 PERIOD ANALYSIS 

In Figure 1 we present the overall 7?-band lightcurve 
(nightly means) which shows that the system has been 



steadily fading from its April 2000 outburst with a decay 
rate of 0.003 magday -1 . After 25 Apr 2001 the rate of de- 
cay slowed and on 12 Jun 2001 the system was found at 
i?=18.650±0.007. The same magnitude was found, within 
the errors, on 26 Jun, so we suggest that this is the true 
quiescent magnitude. It is also consistent, within the errors, 
with the R magnitude from the USNO A2.0 catalogue, where 
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Table 2. Log of spectrometric observations 



Night 


Telescope 


Wavelength 


Resolution 


Exposure 


Number 






range (A) 


(A/pix) 


(*) 


of spectra 


2000 Nov 20 


MMT 


AA4200-7500 


1.1 


1440 


5 


2000 Nov 30 


MMT 


AA4200-7500 


1.1 


1440 


9 


2001 Jan 04 


MMT 


AA4200-7500 


2.50 


1440 


6 


2001 Jan 12 


WHT 


AA5820-7320 


1.47 


1200 


7 


2001 Apr 15 


Bok 


AA4100-7400 


2.75 


1400 


14 


2001 Apr 16 


Bok 


AA4100-7400 


2.75 


1400 


15 


2001 Apr 27 


WHT 


AA6300-6700 


0.41 


340-1000 


23 


2001 Apr 28 


WHT 


AA6300-6700 


0.41 


340-1000 


26 



it is quoted as 7?=18.8 with about ~0.25mag accuracy. Su- 
perimposed on the smooth decay, we also see substantial 
night to night variability. 

Close inspection of the individual light curves shows 
the dominant ellipsoidal modulation of the companion star 
and a distortion wave which progressively moves across the 
orbital light curve. The distortion wave produces dramatic 
changes in the symmetry and minima in the light curves 
(see Fig. ^). As the distortion wave is most probably non- 
sinusoidal in nature we employed the PDM algorithm to 
separate and analyze the different periodicities present in 
the data. 

We detrended the long-term variations by subtracting 
the nightly means from the individual lightcurves. This re- 
moved the low frequency peaks from our PDM spectrum in 
the top panel of Fig. pj (computed in the frequency range 
3 to 30 cycles d _1 with a resolution of 4 x 10 -3 cycles d _1 
and with 25 phase bins). The deepest minima are found at 
5.885 cycles d _1 and 11.769 cycles d _1 which correspond to 
the orbital period (i.e. P =0.169936 d) and the first har- 
monic respectively. A two-component Fourier series (with 
P2 = 2 X Pi; a first-order approximation to the ellipsoidal 
modulation of the secondary star) was fitted to the de- 
trended lightcurve. The fit yielded P orb = Pi = 0.169937(1) 
d and T = 2452022.5122(4), where T corresponds to infe- 
rior conjunction of the companion star. 

The bottom panel in Fig. ^ shows the PDM spectrum 
after subtracting these two Fourier components (P or b and 
Port>/2). A strong signal is found at / = 5.865 cycles d _1 and 
its harmonics. The fundamental frequency corresponds to 
0.17049(1) d, 0.3% longer than the orbital period. Note that 
this perio d is comparable to p eriodicities reported during 
outburst (Uemura et al. 2000b) and hence we interpret it as 
the superhump period caused by a precessing eccentric disc. 
We note this is the first secure evidence of superhumps near 
quiescence in an SXT. 

The light curve of the superhump was obtained through 
the following procedure: first we detrended the overall el- 
lipsoidal lightcurve by subtracting the nightly mean mag- 
nitudes and then phase-folded it into 50 phase bins. This 
was fitted with a simple ellipsoidal model, fixing i — 75° 
and q — 0.04 (see top panel in Fig. ^j. The ellipsoidal fit 
was subtracted from the overall (detrended) lightcurve and 
the residuals were then folded into 50 bins using the su- 



perhump period of 0.17049d (bottom panel of Fig. |3J). The 
lightcurve is clearly non-sinusoidal, and shows a peculiar 
modulation which explains why the PDM spectrum con- 
tained substantial power at the first and third harmonics. 
The shape is clearly different from the superhumps detected 
in other SXTs during outburst (e.g. O'Donoghue & Charles 
1996). However, we note that our data were taken when 
Jl 118+480 was near true quiescence when such features are 
not usually visible. A clear evolution of the superhump from 
single humped modulation to a more complex shape had al- 
ready been observed when the system started its decline^]. 
Peculiar superhump modulations, showing several peaks and 
dips are often detected in cataclysmic variables (CVs) and 
their origin remains unsolved (e.g. V603 Aql, Patterson et 
al. 1993; V503 Cygni, Harvey et al. 1995; H 0551-819, Pat- 
terson 1995). 

Fig. ^ contains nightly lightcurves with superimposed 
model fits of a simple ellipsoidal lightcurve (dashed line) 
and a combination of ellipsoidal modulation plus superhump 
wave (continuous line) . The superhump wave was simu- 
lated by a four term Fourier series fixing the superhump pe- 
riod at P s h=0. 17049 d. The second model is clearly a much 
more accurate representation of the observed lightcurves, 
with x?=l-71 vs x?=3-08. If the superhump modulation is 
due to tidal interactions in an elliptical precessing disc, and 
if the disc's radius shrank as the system faded (as is shown 
in Fig. [j]), we expect the shape and amplitude of the su- 
perhump modulation to change and even vanish if the disc 
radius decreases below the stability radius. 



4 THE ELLIPSOIDAL MODULATION 

The ellipsoidal light curve of the secondary star was phase 
folded and binned into 30 bins using our updated ephemeris 
(see section [j). It can be seen from Fig. ^|that the light curve 
of the superhump modulation is far from simple. It is clear 
that this will distort the ellipsoidal light curve, rendering it 
difficult to interpret since the superhump modulation could 
introduce false features in the light curve. This makes it very 
difficult to interpret any parameters derived from fitting the 
contaminated light curve. However, we can use the observed 



» see the lightcurves in 


http: / /www.kusastro.kvoto- 


u. ac.jp/vsnct/Xray/xtej 11 18-camp.html| 
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amplitude of the modulation and an estimate for the veiling 
in order to to place limits on the binary inclination. 

The veiling for the photometric observations taken in 
April 2001 can be estimated by extrapolating the spectro- 
scopic veiling observed in Jan 2001 (Wagner et al. 2001). 
J1118+480 decreased by 0.30 mags from Jan 2001 to April 
2001. If we assume that this decrease in flux is solely due 
to the accretion disc light fading, then given the observed 
veiling in Jan 2001 of 67%, we estimate the veiling in April 
2001 to be 53%. This value agrees very well with the veil- 
ing estimated from the high resolution WHT spectra taken 
in April 2001; 47±7%, obtained using the standard optimal 
subtraction technique (Marsh, Robinson & Wood 1994) for 
K5-M4 type template stars. 

The large amplitude of the observed modulation 
strongly suggests that the inclination angle is high, but 
the lack of eclipse features at phase 0.0 requires i<82 de- 
grees (also note that no X-ray eclipses were seen during out- 
burst). The difference between the two minima at phase 0.0 
and phase 0.5 in the observed light curve is 0.076 mags. 
However, we have already noted that the observed light 
is heavily veiled (47±7%) and so the true amplitude is in 
the range 0.127-0.165 mags. We have computed the sec- 
ondary star's ellipsoidal modulation assuming a K7v sec- 
ondary star ( r eff =4250 K and logp=5.0; Wagner et al. 2001) 
and q=0.037 ( |Orosz 2001 ) using an irradiated X-ray binary 
model (for details see Shahbaz, 2002). By comparing the 
corrected amplitude with the calculated, we estimate an in- 
clination angle in the range 71-84 degrees. Combined with 
the lack of eclipses this suggests that i lies in the range 71-82 
degrees. 

For more accurate results we would need to model the 
light curves free from the superhump modulation. However, 
deconvolving the two modulations is very difficult because 
in order to do so we have to be confident we have determined 
the 'clean' ellipsoidal light curve. This is only possible once 
we are sure that the system is completely in quiescence at 
which time we expect the amplitude of the superhump to be 
reduced. 



5 H a VARIABILITY 

Our spectra are dominated by strong double-peaked H Q 
emission which exhibits significant time variability both in 
velocity and width. Some nights, for instance, the Ha cen- 
troid and FWHM are modulated with the orbital period 
whereas other nights are consistent with a double sine- wave 
modulation. In addition, the FWHM and line centroid is 
seen to vary by a few hundred kms -1 from night to night. 
Although near quiescence, the accretion disc still contributes 
significantly to the observed optical emission and we see a 
0.5 mag fading in the course of our campaign. This is well 
depicted by the equivalent width evolution of the Ha emis- 
sion line, which shows a continuous rise of ~0.14Aday _1 
due to the continuum decay (Fig. ^|). We also note that the 
Ha line moves dramatically from night to night as is shown 
in Fig. ^ and Fig. Here a three-component Gaussian fit to 
the nightly averaged profiles was performed, consisting of a 
broad base and two narrow peaks. The line is clearly asym- 
metric with the blue peak persistently stronger than the red 
one. We have computed the position of narrow peaks and 



the line centroids given by our fits to the averaged profiles, 
to measure the real line velocity change and not the flux 
weighted centroids. The line centroids move from night to 
night with variations of ~500km/s amplitude. If the nightly 
drifts of the emission lines are due to an eccentric disc pro- 
jecting different amounts of its area at different precessing 
phases, we would expect the variation is modulated with a 
half of the precession period (~26 days). The line centroids 
are consistent with this period with minimum x 2 (showing 
in Fig. top panel), although we note that our spectro- 
scopic sample is not sufficient to rule out a 52 d period (also 
shown in Fig. H). 



6 DISCUSSION 

This paper presents evidence for an eccentric precessing disc 
in J1118+480 during Dec 2000 - June 2001. Although persis- 
tent superhumps have been detected in several extreme mass 
ratio CVs (i.e. Skillman & Patterson 1993) and there is some 
evidence for a superhu mp in the SXT A0620-00 during qui- 
escence ( ]Haswell 199?j ), this is, to our knowledge, the most 
extensive study of a precessing disc in SXTs near quiescence. 
Our lightcurve shows the fingerprint of a distorting super- 
hump modulation with P s h = 0.17049(1) d superimposed 
on the classical ellipsoidal modulation of the secondary star. 
A superhump modulation with P s h = 0.17078(4) d was al- 
ready detected in Ap ril 2000, when the sy stem was at the 
peak of its outburst (Uemura et al. 2000b). This implies a 

which is slower 



superhump period change of P ~ —10" 
but comparable to typical values in SU UMa stars (Warner 
1985). Since the superhump is the beat frequency between 
the orbital and disc precession frequency, the disk precession 
period is, therefore given by P prcc = (P~* - P^ 1 ) « 52 
d. 

The superhump period contains information on the 
mass ratio q. We can use the relation of Mineshige, Hirose 
& Osaki 1992 (hereafter M92) to estimate the mass ratio: 



AP 



Pert 



4 v / TTg 



3/2 



(1) 



where AP is the period excess and r\ is the disc radius in 
terms of of the critical disc r adius for the 3:1 reson ance. If 
we take P sh = 0.17078(4) d flUemura et al. 2000bj) as the 



true superhump period at outburst peak and the empirical 
value r) ~ 0.8 (see M92), we find q ~ 0.028. A similar result 
is obtained using Patterson's (2001) relation, i.e., 

r Psh — Porb-,-1 r 0.37g 



AP = , 

1 Porb J L (l + g) 1 /2 J 

For q in the range 0.04-0.30 he finds 
AP = 0.216(±0.018) <7 



-1 -2.3 



1. 



(2) 



(3) 



for systems with accurate determinations for q (i.e. 77 = 0.8). 
Extrapolating the above expression to Jl 118+480 we find 
q=0.023±0.002. This implies the most extreme mass ratio 
in this class, although we note that for q < 0.025 the 2:1 
resonance rather than the classical 3:1 resonance would be 
the dominant tidal instability (Whitehurst & King 1991). 

The mass ratio has been determined spectroscopically 
by Orosz (2001) using the rotational broadening of the sec- 
ondary stars photospheric absorption lines; q — 0.037 ± 
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0.007. This is comparable to, but less extreme than, the 
mass ratio determined using the superhump period excess, 
but we note that the former is dominated by our assumption 
for 77, thus refined calibrations at the highest and lowest AP, 
in eq. ^ are needed. In this sense J1118+480 is a key system 
in the study of superhumps, since it allows us to extrapo- 
late the 3:1 resonance instability close to the limit where the 
2:1 resonance should start to dominate. In Fig. |§] we show 
the AP-q values for SXTs (O'Donoghue & Charles 1996 and 
this paper) and CVs (Patterson 2001; Molnar & Kobulnicky 
1992 and references therein), including only those systems 
with reliable q determinations. J1118+480 seems to agree 
well with both theoretical AP-q curves for rj = 0.8. 

The nightly mean H Q profiles are shown to move by up 
to ±250 kms" 1 , an effect which has not been seen before in 
any SXT or CV. This behavior is probably driven by the 
changing asymmetric disc brightness distribution also seen 
in our trailed spectra. If the effect we are seeing is the spec- 
troscopic signature of an eccentric disc in J1118+480 pro- 
cessing on a period of 52 d, we would expect a modulation 
consistent with this period as we can see in Fig. [|. In addi- 
tion, the FWHM of the H a line changes by a few hundred 
km s _1 both with the orbital and precession phase. Since 
the FWHM gives a measurement of the velocity dispersion 
in the disc, the observed variability is probably related to 
changes in the disc's shape, since the eccentric disc presents 
differing aspects of its projected area to the observer on both 
the superhump and precession periods. 
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Figure 1. Long-term decay in the lightcurve of J1118+480 from Dec 2000 — June 2001. Only nightly means are shown. Open circles 
represent I-band observations, where a colour correction of R-I=-0.24 has been applied (based on our simultaneous R and I light curves 
of 9 and 10 April 2001). 
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Figure 2. Top panel: PDM spectrum of the whole photometric data base after detrending long-term variations by subtraction of the 
nightly means. Bottom panel: Same but after subtracting the orbital and ellipsoidal (P or b/2) frequencies. The deepest minima are 
marked, corresponding to the superhump frequency (P s h) an d its multiples 2/ and 4/. 
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Figure 3. Top panel: Detrcndcd lightcurvc folded into 50 phase bins and fit with an ellipsoidal model (for i = 75° and q = 0.04). Bottom 
panel: Light curve of the superhump after subtracting the previous ellipsoidal model from the detrended light curve. 
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Figure 4. Orbital lightcurves for a sample of individual nights and best-fit models. Dashed line: ellipsoidal model. Continuous line: 
ellipsoidal model plus supcrhump wave. 
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Figure 5. From top to bottom: Ho centroids (filled circles) and double gaussians peaks (triangles) together with sinusoidal fits of P=26 
days (dashed line) and P=52 days (solid line); mean equivalent width per night with the best linear fit and FWHM of individual spectra. 
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Figure 6. Averaged Ha profiles and three-component fit for each night normalized by the equivalent width, where the positions of the 
peaks and the centroids have been marked. We have computed the phases for the precessing disc using a precession period of 52 d and 
an arbitrary phase 0.0. 
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Figure 7. Upper panels: Observed trailed spectra. Lower panels: Averaged Ha profiles and three-component fit for each night normalized 
by the equivalent width. 
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Figure 8. Relation between AP=(P s ^-P 07 ,(,)/P or 5 and the mass ratio q=M2/Mi for SU UMa stars (filled circles) and SXTs (open 
circles). The open square marks the dipping-bursting X-ray binary V1405 Aqr. The solid line shows the relation given by Patterson 
(2001) and the dashed line the relation by Osaki (1985) for r;=0.8. The vertical line marks the q<0.025 limit for which the 2:1 resonance 
would be the dominant tidal instability. The systems plotted are: WZ Sge, OY Car, Z Cha, IY UMa, HT Cas, DV UMa, V2051 Oph, 
UU Aqr and V1405 Aqr(see Patterson 2001 and references therein); SW UMa, T Leo, V436 Cen, VW Hyi, WX Hyi and TU Men (see 
Molnar & Kobulnicky 1991 and references therein) and VY Aqr (see Thorstensen & Taylor 1997) 
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